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Abstract 

We show that the exactly solved low-temperature branch of the two-dimensional 0(n) model 



' is equivalent with an 0(n) model with vacancies and a different value of n. We present analytic 

2 ■ results for several universal parameters of the latter model, which is identified as a tricritical point. 

These results apply to the range n < 3/2, and include the exact tricritical point, the conformal 

' anomaly and a number of scaling dimensions, among which the thermal and magnetic exponent, 

o ■ 

' the exponent associated with crossover to ordinary critical behavior, and to tricritical behavior 

. with cubic symmetry. We describe the translation of the tricritical model in a Coulomb gas. The 

a\ _ 

' results are verified numerically by means of transfer-matrix calculations. We use a generalized 
OO . 

' ADE model as an intermediary, and present the expression of the one-point distribution function 



in that language. The analytic calculations are done both for the square and the hexagonal lattice. 
PACS numbers: 05.50.+q, 64.60.Cn, 64.60.Fr, 75.10.Hk 
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I. INTRODUCTION 



The 0{n) model is defined in terms of n-component spins on a lattice that interact in an 



isotropic way. Graph expansion [l] of the partition integral of this model leads to a weighted 
sum of graphs, in which every node is connected to its neighbors by an even numbers of 
bonds. In the resulting partition sum, the spin dimensionality n is only a parameter and 
may be varied continuously. For a special choice of the model on the honeycomb lattice, one 
;hus derives a model of non-intersecting loops [2]. Exact results for the universal parameters 
sl as a function of n were obtained for this 0{n) model for two cases, one of them describing 
the critical point, and the other a special (see [j]) case of the low-temperature phase. These 
results are described in the context of the Coulomb gas in a review by Nienhuis 



These two cases o: 



lattice 0(n) model 



m 



universal behavior were also found along two branches of a square- 



together with two different branches describing the critical behavior 
that occurs when 0(n) and Ising degrees of freedom on the square lattice display a joint 
critical point. But it took a long time before an exact description was also proposed j8| for the 
tricritical 0{n) universal parameters in two dimensions as a function of n. The plausibility of 
this description follows from earlier identifications of the fully packed 0(n) loop model with 
the critical Potts model, and of the critical 0(n) model with the tricritical Potts model j^. 
Therefore it seemed plausible that the tricritical 0{n) model is associated with an even higher 
critical Potts model. Such a model is known in the form of the tri-tricritical Potts model, for 
which some universal parameters are known 

were found to agree accurately with numerical estimates of the conformal anomaly and some 
critical exponents of the tricritical 0{n) model on the honeycomb lattice. 

In this paper we provide a more exact basis for this characterization of the tricritical 
0(n) universal parameters, by means of an exact mapping between an ©(ri) -symmetric spin 
model with vacancies, and an undiluted 0{n') model which was exactly solved in Refs. [g] and 
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This part is reported in Sec. [TTl Relations with other models, and the exact solution are 
summarized in Sec. IIIII We include the mapping on the Coulomb gas and derive exponents, 
including the one that is responsible for the crossover to models with cubic symmetry. This 
part is presented in Sec. IIVI In Sec. |V]we define the transfer matrix and apply it to confirm 
numerically some of the predicted exponents. Results and consequences are discussed in 
SeclVD 



II. MAPPING OF THE SPIN MODEL ON A SOLVABLE LOOP MODEL 



As a candidate system to display 0(n) tricriticality, we choose a generalized version of 
the 0(n) spin model on the square lattice studied by Nienhuis and coworkers The 
n-component spins occupy the midpoints of the lattice edges. They are denoted Sj where the 
index i labels the corresponding edge. The integration measure is normalized by J ds = 1 
and the mean length of the spins by / ds{s- s) = n. The distribution is chosen isotropic, 
i.e., in accordance with the 0(n) symmetry. The model also includes face variables tj = 
or 1, which sit on the faces of the square lattice. The index j labels the corresponding 
face. The value tj = corresponds with a 'vacancy' which has the effect of eliminating the 
interactions involving any of the four adjacent spins. The vacancies thus introduce dilution 
and may thus be expected to lead to a tricritical transition, analogous to that in the Potts 
model. 

We write the partition sum as 
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^spin 
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faces j tj=0 
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vert k 



While the spins live on the edges of the lattice, and the vacancies on the faces, the Boltzmann 
weight factorizes into factors W{k) giving the interaction between all variables incident on 
a vertex. A part of the lattice is shown in Fig. [H The local weight W{k) is defined by 
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W{k) = {l+uY^[{p+ {1 - p)ti+2)titi+iti+3 Si ■ Si+i] 



i=l 



+ tit2t3t4Y^[v Si ■ Si+2 +w{Si ■ Si+i){Si+2 " Si+s)] }\/(ti, ^2, ^3, ^4) (2) 



i=l 



where the indices of the s and t variables are defined modulo 4. They describe the position 
of the variables in the local configuration as specified in Fig. [T] around the vertex k. The 
vacancy weight V^(ti, ^2, i^s, ^4) per vertex is specified by 

V{tl,t2, h, ti) = 5t^+t2+ts+t4,A + '^I'^tl+t2+t3+t4,3 + 'y2<5ti+t2+i3+t4,2'^tl,t3 

l+t2+t3+t4,0 (3) 

so that the index of Vi indicates the number of vacancies. The expansion in loop diagrams 
proceeds analogous to Ref. [t], but as a consequence of the added t variables, the loops are 
restricted to the edges that are not adjacent to a vacancy. Furthermore, the term p(l — ^44.2) 
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FIG. 1: Part of the square lattice with spin variables (•), and face variables (o) representing the 
vacancies. The vertices are shown as +. The figure includes the labeling of the variables used in 
the definition of the local weight W{k), when applied to the central vertex in this figure. 



leads to an additional potential for a loop segment that, with respect to a vertex, is diagonally 
opposite to a vacancy. The loop expansion transforms the partition function into 
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^loop 
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(4) 



where the second sum is on all configurations C of closed loops, covering zero or more edges 
of the square lattice, while avoiding edges adjacent to a vacancy. Every vertex is of one of 
ten types shown in Fig. [2l The total number of vertices of type i is denoted iVj, and the 
total number of loops as A^l- The vertex weights Wi are given in Fig. O in terms of the 
parameters that already appear in the spin representation of Eq. ([2]). 



III. EXACT ANALYSIS 



A. Equivalence with the dense loop model 

Consider the exactly solved low-temperature branch of the 0(n) loop model on the square 
lattice, named branch 2 in I?!] and Its partition sum, although of the form of Eq. (jlj), 
is denoted .^dense loop referring to the relatively dense filling of the lattice with loops. The 
weights can be parametrized in terms of the angle 6 G [0,7r/2]. Only the vertices labeled 1 
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FIG. 2: Vertex weights of the 0(n) model with vacancies. The sohd lines represent loop segments, 
and the open circles the vacancies on the faces of the lattice. The presence of a vacancy implies 
that the four surrounding edges are not visited by a loop. The absence of a circle corresponds 
with an occupied face, whose edges may or may not be visited by a loop. The same weights apply 
to rotated versions of the vertices shown here. The spin variables, which sit on the middle of the 
edges, are absent in the loop representation. 

to 4 in Fig. [2] have non-zero weight; their weights are specified as 

W2 = u = h{e) sm{e) 

W3 = v = h{e) sin(3^/4) (5) 
= w = h{e) sm{e/A) 

with h{9) = l/[2 sin(6') cos(36'/4) + sin(36'/4)]. The weight of the loops (or the dimensionahty 
of the spins) is 

n' = -2cos(20) . (6) 

We use n' because we wish to reserve n for another choice, in which all weights Wk are 
non-zero. 
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FIG. 3: The adjacency diagram A having q branches. For g = 0, 1 and 2, it reduces to the diagrams 
A2, A4 and Eg respectively. 

B. ADE models 



Here we construct an alternative representation of the loop model partition sum. Follow- 
ing Pasquier Q , the loops are interpreted as domain walls in a configuration of discrete 
variables living on the faces of the lattice. These variables take values, corresponding to the 
nodes of a graph A called the adjacency diagram. In this paper we consider the family of 
graphs shown in Fig. [31 but the discussion in this section is general and the figure can be 
seen as an example. We call this model an ADE model after the classification of adjacency 
diagrams. Neighboring faces not separated by a loop carry the same value. If they are 
separated by a loop, their values are adjacent in A (hence its name). 

Associated to A we introduce the adjacency matrix A with elements Aj, where i and 
j represent nodes of A. The elements are defined as Aij = 1 if z and j are adjacent, and 
Aij = otherwise. Of the 2(1 + q) eigenvectors of A, four are symmetric under permutation 
of the q branches. In the symmetric subspace the eigenvector equation can be written as 
that for the right hand eigenvector of the 4x4 matrix Aij 
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(7) 



The first two elements of the eigenvectors of A correspond to the nodes a and b respectively. 



and the last two with the nodes of type c and d respectively (see Fig. [3]). We choose 
the number of branches of the diagram ^ to be g = {n' — l/ra')^, so that the symmetric 
eigenvalues of A are given by 

Ao = n' 
Ai = l/n' 

A2 = -l/n' (8) 
A3 = -n', 

and for each of these eigenvalues, the elements of the corresponding eigenvectors are 

(^^ S^,) = (A, - 1/A, , A,(A, - 1/A,) , A, , 1) . (9) 



Besides these there are 2(g — 1) eigenvectors 
branches of A. They have eigenvalues 

A4 = 
As = 



antisymmetric for interchange of two of the 

1 (10) 
-1, 



each (g — l)-fold degenerate. 

Of the eigenvectors of A, generally denoted as Sj, we omit the upper index for the case 
/i = 0. For n' > 1 this is the Perron- Frobenius eigenvector, with eigenvalue Aq = n' > A^ 
for /i = 1, 2, 3. 

We write the weight of the corresponding ADE model in terms of a product of local 
weight factors as 

(C \ 7bcnd/27r 
^] ■ (11) 

There is a factor W{k) for each vertex depending only on the local configuration of domain 
walls, and a factor for each turn of the domain wall which also depends on the states of the 
faces on the inside (i) and outside (j) of the loop. The bending angles 7bend are counted 
positive where the loop bends inwards, so that the sum of the bending angles along a loop 
is +27r. Thus the weight of an entire closed domain wall is 

Wioop = A^,^, (12) 



and still depends on the state i inside and j outside the domain wall. When for a fixed 
configuration of domain walls the sum over compatible state configurations is performed, 
each closed domain wall thus contributes a factor equal to the largest eigenvalue of A, i.e. 
n', just as in the 0{n') loop model. This confirms that the partition function of the ADE 
model satisfies 

ZaDE = X! W^ADE = -^dense loop , (13) 

g {s}\g 

where the first sum is on all loop configurations and the second one on the configurations 
of ADE variables compatible with the loops. The suffix of the right hand side emphasizes 
that the loop model is in the dense phase, and does not permit vacancies, which will be 
introduced later. Since in Z^ense loop the variable n' and therefore q only enter as parameters 
we may vary them continuously, thus representing the continuous n'-weight loop model. 

When the models live on a torus rather than in the plane, there may be loops that wind 
the torus. In the loop model these typically have the same weight n' as the contractible 
loops. However, in the ADE model, the corresponding domain walls have a net bending 
angle equal to zero. This implies that they carry the weight Aj^k rather than Aj^kSj/Sk- The 
summation over the states of the domains then reduces to taking the trace of a power of 
A, equal to the number of non-contractible loops. The result is that all the winding loops 
have the same weight, equal to an eigenvalue of A, which should then be summed over these 
eigenvalues. We conclude that the ADE model on a torus corresponds to a loop model in 
which the winding loops receive special treatment. Or alternatively the loop model partition 
sum with all loops weighted equally is the largest sector of the ADE transfer matrix. 

C. Correlation functions 

We will now calculate the one-point distribution (IPD) P{k) of the ADE model, i.e., the 
probability that a face is in state k. Consider a loop well inside a large lattice. We assume 
that the IPD is unaffected by the presence of the loop (or any other loop). In other words we 
assume that the IPD conditional on the presence of a loop is the same as the unconditional 
IPD. That this is plausible follows from the calculation of the partition sum above: the 
contribution to the partition sum of a particular domain is independent of the domains it 
is contained in, and it is independent of all the domains it contains, once the state of these 
domains has been summed over. 

o 



The conditional probability P{k\j) that the inside domain of a loop is in state k, provided 
the outside domain is in a given state j, is determined by Eq. f|T2|) as P{k\j) = AkjSk/{n' Sj). 
Thus we find the joint probability P{k,j) that the outside of a loop is in state j and its 
inside in state k as 

p{k,j) = p{mk\j) = p{j)Ak, 4k- ■ (14) 

n 

Summation on j now yields the probability that the inside domain is in state k, which should 
be equal to P{k): 

Y^P{k,j) = P{k). (15) 
i 

Using the symmetry of A, one finds the unique (normalized) solution to this consistency 
condition as 

P{k) = — ^ . (16) 

An approach alternative to the condition that the loop considered is well inside a large 
lattice, is to consider a bounded lattice of arbitrary size, with the faces on the boundary 
all in the same state, with Eq. (fT6l) as the probability distribution for that state, the ideal 
fixed boundary condition. Then by induction the same distribution holds for the domains 
separated from the boundary by one domain wall, and so on recursively to the innermost 
domains. It is then assumed that in the thermodynamic limit the boundary condition should 
not matter, well away from the boundary. 

Consider the function S^/Sk, i.e., the ratio of an arbitrary eigenvector S'^ and the Perron- 
Frobenius eigenvector 5*. If this function is part of a correlation function (■ ■ • Sji/Sk ■ ■ ■), 
where k is the state of a given face, it effectively changes the weight of the loops surrounding 
the face. This is easily seen in the expression ( IT2|) : the factor S^/Sk replaces the numerator 
by S*^, so that the weight of the loop becomes that of the corresponding eigenvalue, A^, 
as long as they do not surround other operator insertions. We will call these functions 
weight-changing operators. 

A more interesting result |l^ comes from the two-point function 

(f f > • 

j and k being the state of two arbitrary faces. The weights of the loops surrounding either 
of these faces but not the other is changed into the respective eigenvalues and A,^. 
corresponding to the eigenvectors S'^ and S'^. Now consider the innermost domain wall that 



surrounds both faces. After the states of the domains nested inside it are summed over, the 
weight governing the state of the final domain is 

, (18) 

where k is the state of the surrounding domain. This can be expanded as a hnear combination 
of all eigenvectors: 

s^s^ = s,Y.c;.s^, (19) 

where, provided the eigenvectors are normalized, 

cA* qiy OK 

c;. = T.^^- (20) 

Apparently the combination of two operators labeled fi and u look from a distance like a 
linear combination of operators k. 

These structure constants of the operator product expansion, or fusion rules, may readily 
be calculated explicitly for the diagrams in Fig. [3l but here we only note that they are 
symmetric in /i, u and k, and that they vanish if one of the indices corresponds with the 
largest eigenvalue, and the others two differ. This implies that the two-point correlation 
function of two different weight-changing operators vanishes in the thermodynamic limit. 
Obviously these fusion rules may be used just as well in correlation functions of more than 
two operators. 



D. Equivalence with the 0(n) model with vacancies 

We interpret the extrema of A as vacancies, that is the nodes a and those of type d. 
Thus, there may be g + 1 types of vacancies in an ADE configuration. However, in any 
configuration the type of each vacancy is fully determined by the neighboring domains. 
Thus, it is sufficient to specify the g + 1 states for the non- vacant domains in order to fully 
describe an the ADE configuration (with the exception of the completely vacant state). 

Now we identify the domain walls between domains in state b and in states of type c as 
loops, so that we have a loop model with vacant faces. Note that for any given configuration 
of vacancies this loop model is much like that described in Sec. IIII Bl but on a restricted 
lattice, from which the vacant faces are omitted, and with a reduced adjacency diagram, of 
only the full nodes in Fig. [31 Note that the eigenvectors of the reduced adjacency matrix, up 

1 n 



to normalization, are the same as the eigenvectors of the total adjacency matrix, restricted 
to the nodes b and those of type c. The eigenvalues are ^^/q■ Following the arguments used 
before, we obtain a loop model with loop weight n = y/q, with weights that follow from the 
original ADE model with the complete adjacency diagram and the entire lattice. Thus the 
successive transformations are 



^ dense loop 



•^loop+vac. ) (21) 



in which the right hand side is the partition sum of a model with loops and vacancies. The 
weights of this loop model are, with reference to Fig. [21 given by 
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(22) 

respectively. The two terms in IVg arise from the two orientations of the type-4 vertex of 
the O(n') loop model. 

We note that the weights are completely given by the configuration of loops and vacancies, 
irrespective of the type of vacancy. Furthermore, any configuration of loops and vacancies 
consisting of the local vertices in Fig. [2]is possible for the adjacency diagram in Fig. [31 These 
properties are not generic for any adjacency diagram, and are the basis of our choice of the 
diagram in Fig. [3l together with the fact that is contains a continuously variable parameter 
controlling the eigenvalues of its adjacency matrix. 
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FIG. 4: Vertex weights of the honeycomb 0(n) model with vacancies on the faces. Sohd Unes 
represent loop segments, and circles the vacancies on the elementary hexagons. A vacancy excludes 
the six surrounding edges to be visited by a loop. Rotated versions of these vertices have the same 
weights. The spin variables, which sit on the vertices, play no part in the loop representation. 

E. The 0(n) model on the honeycomb lattice 

A similar 0(n) spin model with vacancies on the faces can be defined on the honeycomb 

I I 

lattice, see e.g. Ref. '8'. Here we also include interactions between the vacancies, described 
by the three vertex weights fi, and where the indices show the number of vacancies 
adjacent to the vertex. 

The transformation into a loop model partition sum proceeds the same as for the square 
lattice, and leads to the form of Eq. (jlj) but with only five independent vertices. They are 
shown in Fig. IU together with their weights. For the simplified case w i = f 2 = = 
without vacancies, and with the special choice 
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1/V2 ± V2 
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(23) 
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jj. The 



(which is different from that for the square lattice), this model is solvable 
high loop density branch of this model corresponds with the minus sign. 

The mapping of the low-temperature 0{n') model on the dilute 0(n) model can be 
performed analogously for the honeycomb lattice, and leads to the following vertex weights 
for the model with vacancies: 



1 



u 



The adjacency matrix, and thus its eigenvector components, are the same as for the square 
lattice model, Eq. ([9]). Also the loop weight 

n = n' - — , (24) 
n' 

is the same. 

F. Interpretation 

It remains to be shown that the constructed model of loops and vacancies is a tricritical 
0(?7.) model. To this purpose we obtain its conformal anomaly c via the equivalence with 



the low-temperature 0{n) model of Eq. ([5]), for which 



m 



6(1 -oY 1 

c = 1 - ^ — , 2 cosfTT^) = -n', 0<^<1, n = n' . (25) 

9 n' 

This combination of c and n does not agree with the known critical and low-temperature 

n 

0(n) universality classes [3|. Since it is known that the introduction of vacancies can lead 
to tricriticality, this already suggests that the model defined by Eqs. ([1]), (1221) and is 
tricritical. Further justification will be given below. 

IV. UNIVERSAL PROPERTIES AND MAPPING ON THE COULOMB GAS 
A. The conformal anomaly 

Eq. (125|) shows that the model with vacancies derived above does not fit the critical and 
low-temperature 0{n) universality classes It does, however, precisely match the result 
for c inferred in Ref. [sl for the tricritical 0(n) model, namely 

6 TT 1 

c = 1 , 2 cos — — = A, m>l, n = A - — . (26) 

m(m + 1) m + 1 A 

The parametrizations of n in Eqs. ( !25l) and ( !26l) imply that n' = A. This provides a 
confirmation that the 0(n) loop model defined by the vertex weights of Eq. ( l22l) and the 
loop weight of Eq. ( |2^ is a tricritical 0(n) model. For the Ising model n = 1 and the self- 
avoiding walk n = this proposal reproduces the the known tricritical values of c = 7/10 
and c = respectively. 



B. Relation with the Coulomb gas 



The Coulomb gas offers a powerful tool to calculate critical exponents, provided the 
Coulomb gas coupling constant g is known for the model under investigation. The coupling 
constant is related [isl, [3] to the conformal anomaly c according to Eq. fl25l) . In combination 
with Eq. ( 126|) this yields as a function of m which is two- valued. For the tricritical 0(n) 
model one has g = m/ (m+1), in agreement with the conformal classification of the tricritical 
Ising {n = 1) model. In the Coulomb gas language, the conformal anomaly is represented 
by means of a pair of electric charges eo as follows: 



1-^, eo = l-g. 



n 



—2 cos Tig . 



Combination with another pair of charges ±e^ yields a scaling dimension 



(27) 



(28) 



n 



which reproduces the dimensions listed in Refs. 18| andl9| for the tricritical 0(n) universality 
class for charges according to 



cos Tie, 



(29) 



The An follow from the eigenvalues of the adjacency matrix |12l . Il3l ]. which is the same as 



M]- Six eigenvalues were already listed in Ref, 



that used for the tri-tricritical Potts model 
and in Eqs. ([8]) and ffTOl) . The leading scaling dimensions follow as 
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with 
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These exponents are associated with the weight-changing operators Sjl/Sk, discussed in 
Sec. nil CI This implies, see Eq. iQ, that the X2 is associated with the enhancement (or 



suppression) of vacancies: its eigenvector has a different signs in the vacant and non-vacant 
states. It follows that X2 is one of the thermal exponents, leading or otherwise. The 
exponent Xi belongs to the operator that breaks the pseudo-symmetry between the left and 
right hand side of A. It changes the weight of the loops from n to —n, but is not meaningful 
in the local 0(n) spin version of the loop model. One can arrange this effect, however, by an 
operator that terminates a seam across which the interactions have the opposite sign. It will 
appear later as the interface exponent. The exponent X3 combines the effect of the prior two 
operators. Finally the the operator associated with the exponent X4 breaks the symmetry 
between the q branches, while that of X5 combines this with the enhancement of vacancies. 
Again, these operators only show up in models where the nodes of A are meaningful. In our 
0{n) models, the diagram is only used to give the proper weights to various configurations. 

In addition to the dimensions describing singularities associated directly with the ADE 
model, we consider exponents describing crossover phenomena due to a possible perturbation 
of the 0(n) symmetry. Such perturbations can, in the Coulomb gas language, be associated 
with pair of electric charges 1 — g combined with a pair of magnetic charges whose magnitude 
depends on the type of perturbation [sj. A pair of magnetic charges ±k, associated with 
the defect of 2k lines coming together, then corresponds with a scaling dimension 

^^''-^-2^^ 2 ^^^^ 

C. Specific exponents 

The magnetic exponent 

In the 0(n) loop model, the magnetic correlation function is represented by configurations 
containing a single loop segment connecting the correlated points. In the Coulomb gas 
language it corresponds with magnetic charge k = 1/2. Then, Eq. (!32|) yields 

X, = l-1-^. (33) 

This is in agreement with an earlier conjecture 20|] and with numerical results for a related 
model with vacancies on the honeycomb lattice [8]. After the mapping onto the model with 
vacancies, the defects, i.e., the end points of the loop segment, can only sit on non-vacant 
sites, and not in the regions occupied by vacancies. But otherwise, the defects have the same 



physical effect, and therefore the same exponent (!33|) . Therefore the magnetic exponent of 
the tricritical 0{n) model is the same as that of the low-temperature 0{n') model. 
The temperature exponent 

Above we already identified the exponent X2 as a thermal exponent. This agrees with 
identification on the basis of numerical evidence, in 8(] as, in fact, the leading thermal 
exponent. Effects described by this exponent are absent in the thermal properties of the 
low-temperature O(n') model. This fits well in the interpretation [3] that the latter model 
is confined to the critical subspace of a larger parameter space. Such 'unphysical' O(n') 
exponents are however known to reappear in correlations on dilute 0(n) loop configurations 



2l| . We further remark that the thermal exponent of the O(n') model, associated with 
the suppression and enhancement of loops, should be another thermal exponent; its value, 
4/(7 — 2, however indicates that it is irrelevant. Another option to find a second thermal 
exponent, comes from the observation that the value of t in Eq. (13T|) is only the smallest 
solution of the inverse cosine. The next leading exponent is obtained by replacing thy 2 — t 
in X2, so that 

A, X. ^ P-"'-;'-^'^ (34) 

Other exponents 

The introduction of a 'seam', i.e., a row of antiferromagnetic bonds, such that the bonds are 
perpendicular to the row, leads to a change of the partition sum described by an 'interface 
dimension Xm- While in Ref. [g the identification X^ = Xi was made on the basis of 
numerical evidence, here we can make the identification by inspection of the corresponding 
operator. A seam along the length of the cylinder changes the sign of the non-contractible 
loops. That corresponds precisely with Ai = 1/n' in Eq. ([8]) because inversion of n' results 
in a change of sign of n. 

In the 0{n) spin model, other interface exponents can be constructed by the introduction 
of a cut across which the spin sis identified with R-s, where R is an operator in the orthogo- 
nal group 0(n). The weight of the loops crossing this cut is then equal to J ds{s ■ H ■ s). By 
varying R this can take any value. The corresponding exponent is obtained by the relation 
A - 1/A = Jdss-R-s, and Eqs. ([28]) and ([29]). 

A cubic perturbation of the 0{n) symmetry can be represented by magnetic charges 



k = ±2m Eq. ([32]) |3|, [2^ which yields 



a 



AV2=l-^ + |. (35) 

It is the exponent that describes the crossover when a cubic symmetry breaking is intro- 
duced. Also when intersections between the loops are permitted, this exponent governs the 
crossover to another universality class. This indicates the fact that the results in this paper 
are applicable exclusively when intersections are prevented, by the specific choice of the 
Hamiltonian. 



V. NUMERICAL VERIFICATION 



For the construction of the transfer matrix we choose the usual geometry of a model 
wrapped on a cylinder, such that one of the lattice edge directions runs parallel to the axis 
of the cylinder. The transfer-matrix method used here is based on that of Ref. |7|, including 
the sparse-matrix composition. The main modification is the generalization of the set of 
connectivities used in Ref. Q to include the specification of the vacancy variables on the 



faces. 



A. Enumeration of the connectivities 



We consider the model of Eq. (jlj) on a cylinder with a circumference of L lattice units. 
The cylinder is has an open end such that there are L external edges, which may or may 
not be covered by segments of incomplete loops. The connectivity specifies the following 
information: (1) which of the faces at the end of the cylinder carry vacancies; (2) the way 
in which pairs of covered external edges are connected by incomplete loops of C. These 
connectivities are subject to the restriction that the loop segments cannot be adjacent to 
a vacancy. Each connectivity can be fully specified by a row of integers {ii,i2, - ■ ■ , ii) such 
that 

k = im > if and only if edge / is connected to edge m 

ik = if and only if edge k is not visited by a loop segment and the 

< 

face to the right of k is occupied 
in = —1 if and only if face to the right is vacant. 

1 rr 



The positions of the vacancies can simply be coded by means of an L-bit binary number 
{pi,P2, ■ ■ ■ ,Pl) with value P = T.^^^pk^''"^ + 1. For a given /5, we no longer need those ik 
that sit adjacent to a vacancy. After dropping these ik from {ii, ■ ■ ,iL), let (ji, ^2, ■ ■ ■ ,ju) 
denote the remaining sequence of length u. This sequence can be coded by means of an 
integer cr(ji, j2, ■ ■ ■ ,ju) in the range I < cr < ttu- The actual values of a and of are given 
in Ref . Q. Let 



be the number of connectivities whose binary vacancy number is smaller than /3, where 
u{a) is the number of dangling edges which are not adjacent to a vacant face in the face 
configuration a. Then, the integer that codes the connectivity with vacancies is 



A decoding algorithm, that constructs a sequence ■ ■ ■ , ^l) given the integer 7, was 

constructed using similar methods. 

B. Numerical calculations 

Several eigenvalues of translationally invariant (zero-momentum) eigenstates of the trans- 
fer matrix were computed for a limited range of system sizes L < 16, as follows: 

1. The largest eigenvalue A^^"* in the 'even sector', which means that the transfer matrix 
operates in the space of connectivities whose dangling bonds occur only in connected 
pairs. 

2. The second largest eigenvalue A^^-* in the same sector. 

(2) 

3. The largest eigenvalue A^^ in the 'odd sector', which means that the transfer ma- 
trix operates in the space of connectivities with, apart from dangling pairs of bonds, 
precisely one dangling bond that is single. 

4. The largest eigenvalue A^^ in the even sector of the transfer matrix of a model with 
a 'seam'. The seam modifies one row of bonds. These bonds are perpendicular to the 
axis, while the row itself is parallel to the axis. All edges of this seam contribute a factor 
— 1 to the Boltzmann weight, if covered by a loop segment. In actual calculations, this 




7(^1, ^2, ■ ■ ■ , «l) = MP) + (^Ul,j2, ■■■,ju)- 



(36) 



1 Q 



is realized by changing the sign of some of the vertex weights of Fig. [2] and Eq. (!22l) . 
for those vertices that are immediately to the left of the seam. 

The finite-size data for the largest eigenvalue A?'* determine the free energy density, from 



which we estimated the conformal anomaly c [18|, ll9|. The ratio A^^^/A^^'' defines the corre- 



ation length of the energy-energy correlation function. Using Cardy's conformal mapping 



231] of an infinite cylinder on the infinite plane, one can thus estimate the temperature di- 
mension Xf. Similarly, A^^^/A^^"* is used to find the magnetic dimension X^- Finally the 
ratio A^'*/A^°'' yields the so called interface exponent Xjnt- All of the quantities c, Xt, Xh, 
and Xint were already described exactly as a function of n, and verified numerically, see 
Ref. ^ and references therein. The present numerical analysis is aimed at confirming that 
the present model describes the tricritical 0(n) model. The numerical analysis follows ba- 
sically the hues of Refs. Q and 1^; see also b^. The final estimates are hsted in Tabs. [T] and 
Hn They agree convincingly with the analytic expressions listed in Sec. IIVI whose values are 
also included in the tables. 



VI. DISCUSSION 

The present tricritical OJn) model appears to belong to the same universality class as a 
loop model defined in Ref. M- The latter model was defined as the surrounding loop model 
of the critical g-state random cluster model on the square lattice. It is possible to apply 
the same method as used above, namely to use the ADE interpretation and to restore the 
loops except those surrounding the vacancies of type a and d, to the latter loop model. We 
have chosen the present formulation, based on the low-temperature 0(n) model of branch 

n p 

2 defined in Refs. |6| and |7|. This is more natural in the sense that it allows for sites that 
are neither visited by a loop, nor adjacent to a vacancy. The relations between the various 
models, as constructed and listed in Sees. [TTl and llllt are summarized by 

Zsp'm -^dense loop -^ADE ^ -Z^loop+vac. ^ -^dilute spin • (^'^) 

The last step follows from the general equivalence formulated in Sec. [TTl 

As found in Ref. sl, the introduction of vacant faces in the honeycomb 0{n) model leads 
to tricriticality when the fugacity of the vacancies is sufficiently large. No vacancy-vacancy 
couplings were introduced. Numerical work on the square lattice 0{n) model for n = 1 



TABLE I: Conformal anomaly c and interface critical dimension as determined from the 
transfer- matrix calculations described in the text. Estimated error margins in the last decimal 
place are given in parentheses. The numerical results are indicated by '(num)'. For comparison, 
we include theoretical values. 



n 


c (num) 


c (exact) 


Xm (num) 


Xm (exact) 


-2.0 


-0.99155 (1) 


-0.9915599 


-0.20179901 (1) 


-0.201799000 


-1.75 


-0.91099 (1) 


-0.9109986 


-0.17697229 (2) 


-0.176972272 


-1.50 


-0.81973 (1) 


-0.8197365 


-0.15164470 (2) 


-0.151644706 


-1.25 


-0.71646 (1) 


-0.7164556 


-0.1259301 (1) 


-0.125930086 


-1.00 


-0.59999 (1) 


-6/10 


-0.100000000 (1) 


-1/10 


-0.75 


-0.46962 (1) 


-0.4696195 


-0.07409548 (2) 


-0.074095457 


-0.50 


-0.32528 (1) 


-0.3252829 


-0.048531921 (1) 


-0.048531921 


-0.25 


-0.16799 (1) 


-0.1679953 


-0.023691688 (1) 


-0.023691689 

















0.25 


0.175264 (1) 


0.1752630 


0.02211104 (1) 


0.0221110351 


0.50 


0.353480 (1) 


0.3534792 


0.042235700 (1) 


0.0422356998 


0.75 


0.529949 (1) 


0.5299489 


0.060000362 (1) 


0.0600003616 


1.00 


0.700000 (1) 


7/10 


0.07500000 (1) 


3/40 


1.25 


0.85897 (1) 


0.8589769 


0.086505216 (2) 


0.0865052157 


1.50 


1.00000 (1) 


1 


0.08801923 (5) 


0.0880192310 



revealed a peculiar difference with the honeycomb 0{n) model. No tricritical point was 
found when vacancy- vacancy couplings are absent. Instead, a multicritical point resembling 
that of branch 3 of Ref. [71 was found. The physical interpretation of this multicritical point 
is that the 0(?7,) critical line merges with an Ising critical line, where auxiliary variables in 
the form of dual Ising spins undergo a phase transition. A qualitative difference with the 
model described by the vertex weights of Eq. (l22i) is that the vacancies attract each other 
in the latter model. 

A comparison of the numerical results for the present model with those for the tricritical 
honeycomb model 0{n) studied in Ref. 18! shows a conspicuous difference in the estimated 



TABLE II: Temperature critical dimension Xt and magnetic dimension X^ as determined from 
the transfer-matrix calculations described in the text. Estimated error margins in the last decimal 
place are given in parentheses. The numerical results are indicated by '(num)'. For comparison, 
we include theoretical values. 



n 


Xt (num) 


Xt (exact) 


Xh (num) 


Xh (exact) 


-2.0 






-0.0951628 (1) 


-0.0951627339 


-1.50 




0.709784688 


-0.0876432 (1) 


-0.0876431495 


-1.25 


0.4814737 (2) 


0.481473928 


-0.0790909 (1) 


-0.0790908776 


-1.00 


0.39999999 (1) 


2/5 


-0.05833333 (1) 


-7/120 


-0.75 


0.3446680 (2) 


0.344668096 


-0.045889544 (1) 


-0.0458895426 


-0.50 


0.3039307 (2) 


0.303930873 


-0.031982842 (2) 


-0.0319828413 


0.00 


0.2500000 (1) 


1/4 








0.25 


0.2324956 (1) 


0.232495729 


0.017729518 (1) 


0.0177295181 


0.50 


0.2192386 (1) 


0.219238626 


0.03627658 (1) 


0.0362765827 


0.75 


0.2088742 (1) 


0.208874121 


0.05539746 (1) 


0.0553974632 


1.00 


0.20000001 (1) 


1/5 


0.07500000 (1) 


3/40 


1.25 


0.19068002 (1) 


0.190680043 


0.09549715 (1) 


0.0954971419 


1.50 


0.16844985 (5) 


0.168449854 


0.125000000 (1) 


1/8 



accuracies. This difference can be explained from tlie way in wliicli ttie two different sets of 
tricritical points were found. For tlie lioneycomb lattice 0(n) model of Ref. Is, tlie tricritical 
points were determined numerically in a small parameter space. From the perspective of the 
renormalization theory, this procedure yields a rather arbitrary tricritical point in the sense 
that the irrelevant fields are non-zero in general, and thus introduce corrections to scaling. In 
contrast, the exact equivalence of the present tricritical square-lattice 0(n) model with the 
0(n) low-temperature branch indicates that the leading irrelevant field vanishes, since the 
equivalent 0{n) low-temperature branch is characterized by the vanishing of its irrelevant 
temperature field. As a result the corrections to scaling are suppressed, and the apparent 
finite-size convergence improves drastically. 

Since it is widely believed that the universal parameters describing the critical state are 



determined by the symmetry of the model, the dimensionahty, and the range of interaction, 
it seems plausible that the tricritical model presented above serves as a representative of 
the generic 0(n) universality in two dimensions. Indeed the spin-spin interactions defined 
in Sec. [IT] contain only scalar products, which satisfy the 0{n) symmetry. 

However, in this case the 0(n) symmetry of the spin model, is not a secure guide for the 
universality class. This is because, like in the dense loop phase of the pure 0{n) model, 
intersections are relevant . The same applies to the tricritical point reported here: the 
exponent associated with crossing loops is the same as that of cubic symmetry breaking. 
Recently, Jacobsen et al. proposed that the low-temperature phase of the generic 0(n) 
model is described by the intersecting loop model proposed in [25| and since called Brauer 



model 



26|. 



It is interesting to note t 



the critical 'branch 



i-aa 



lat the mappings described in Sec. [TTT] can also be applied to 



ll| of the square- lattice 0(n) model. Just as branch 1 is the 



analytic continuation of branch 2, we can continue the tricritical branch through the 'end 
point' n = 3/2, g = 1 to g > 1. The weights for branch 1 are also given by Eq. (JSj), but 
instead with 7t/2 < 6 < n. The relation between n and n' remains the same, but the vertex 
weights as specified by Eq. change, and the relation between the Coulomb gas coupling 
and the conformal classification parameter m is no longer g = m/(m + 1) but becomes 
g = (m + l)/m > 1, while it relates to n' as n' = 2cos(7r/m) (see e.g. Ref. Q). For = 1 or 
n' = (1 + V5)/2 one thus finds a higher critical Ising model that is to be compared with the 
m = 5 model in the series of Andrews et al. [27!|. For n = or n' = 1 the model displays 
an Ising-like critical point. The resulting branch of multicritical points can thus be seen as 
a generalization of the m = 5 Ising-like model for n = 1 to continuous values of n, i.e. the 
point where the tricritical point itself turns first-order. 
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